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Abstract 

Recent lattice computations showed how the approach dubbed "refined Gribov-Zwanziger sce- 
nario" is in very good agreement with data and the gluon propagator fits them very weh. This 
propagator can be described as a finite sum of Yukawa propagators and can be obtained introducing 
some condensates representing the contribution of the vacuum of Yang-Mills theory. The values of 
the condensates are arbitrary and are obtained through lattice data. This kind of structure of the 
propagator is in agreement with the one we obtained using a different analytical approach but with 
the substantial difference that here all the physical parameters are properly fixed. We show that 
our approach can properly fix all the values of the condensates in the refined Gribov-Zwanziger 
scenario giving a complete validation to both techniques. This will provide an interesting view to 
the gluon propagator and a set of values to be experimentally determined. 
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I. INTRODUCTION 



In recent years, lattice studies of the two-point functions made clear their structure l|-l3|. 
These results have shown how a previous view 4, |5| about gluon and ghost propagators 
for a pure Yang-Mills theory was diverging in giving an explanation to their behavior. The 
question that was put forward by Gribov and then improved by Zwanziger [3| showed 
that the gluon propagator should go to zero as momenta go to zero and the ghost propagator 
should go to infinity faster than a free propagator. These conditions granted that positivity 
was maximally violated and the theory is confining. The new paradigm emerging from 
lattice computations showed instead a rather different situation: Gluon propagator reached 
a finite non-null value at zero momenta and the ghost propagator was indeed behaving as 
that of a free particle. This situation is in good agreement with the idea of a massive gluon 
(a mass gap) for Yang-Mills theory but appeared to conflict with the initial ideas about 
confinement. 

From a theoretical standpoint, the pioneering work of Cornwall in the '80 was claiming 
for the idea of a non-null propagator at zero momenta js]. At the start of this century this 



idea was revamping jo], |10| wit 
theoretical work that followed 



1 numerical analysis of Dyson- Schwinger equations and the 



ll| to support this conclusion. 



Between different theoretical approaches that are emerging in order to give an under- 
standing to the unexpected behavior of the Green functions of Yang-Mills theory, the idea 
that condensates are playing some relevant role in the vacuum of the theory was initially put 



forward 



12j |. This scenario has been recently improved in such a way to completely improve 



Gribov-Zwanziger scenario in a "Refined Gribov-Zwanziger" scenario (RGZ 



mg 



13 



or the follow- 
These 



16 



- |l5| that reached a substantial agreement with lattice computations 
comparisons with lattice data are a proof that condensates, that are generally postulated in 
this approach, indeed exist and this way to analyze the low-energy behavior of Yang-Mills 
theory is sound. 

The structure of the gluon propagator that emerges from these studies has the form of a 
sum of free propagators of massive particles. This idea was firstly put forward by Migdal in 

jy Bochicchio for large-N Yang-Mills theory give strong support 



the 70 



19|. Recent works 



to Migdal's conclusion 



20- 



22|. 



In view of these relevant conclusions, it appears really important to derive from first 
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principles the values of the condensates and make the analysis complete from numerical to a 
theoretical one. In order to accomplish this task, we will use a technique devised by us that 
has the advantage to have all the parameters properly derived from the theory and to ge t 
the proper structure of the gluon propagator as a sum of free massive propagators 23N25|. 
In this way, we will see this as a validation to both these approaches being them each other 
consistent. The idea behind this technique is to take for granted, at some stage, that the 
emerging data from lattice about the running coupling of the theory are just displaying a 



trivial infrared fixed point for a pure Yang-Mills theory [26| . Then, by proving the existence 
of a class of classical solutions to Yang- Mills theory, this behavior is shown to emerge and, 
combined with a Gaussian form of the generating functional in the in the infrared limit, we 
can build up a consistent framework to compare with RGZ scenario. 

The paper is so structured. In secjlT] we discuss the gluon propagator as given in RGZ 
scenario and derive it from instanton solutions obtained from the classical equations of 
motion of the Yang-Mills theory, starting from lattice data showing the running coupling 
to go to zero lowering momenta. In sec JIIII we derive the condensates of RGZ scenario 
in a closed form proving also that the agreement is consistent with a sum of two Yukawa 
propagators as also is emerging in RGZ studies and fits with lattice data. Finally, in sec JIVI 
conclusions are given. 

II. GLUON PROPAGATOR 



A. Refined Gribov-Zwanziger scenario 



a 



The initial formulation of the Gribov-Zwanziger scenario was given in 1989 |7|. A com- 
pletely renormalizable Lagrangian of Yang-Mills theory was formulated that agrees at all 
orders with Gribov prescription to restrict the path integral to the first Gribov region. This 
was achieved with the introduction of a couple of bosonic fields ((/:>^, ip^) and a couple of 
fermionic ones {u^^LiJ^^). This Lagrangian provides the following prescriptions for the gluon 
propagator in the Landau gauge in the infrared limit for SU(N) 

with 7 is a constant named Gribov parameter. This propagator goes to zero lowering 
momenta and so is in disagreement with lattice computations. Similarly, this scenario gets 



a ghost propagator infrared enhanced as it goes to infinity faster than the free propagator in 
the same hmit, again in disagreement with lattice data. But this scenario can be reconciled 
with numerical data when one notes that the effective potential in the Gribov-Zwanziger 



theory permits the presence of condensates 
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can be written down in the Landau gauge as 



-llSll . Then, the most general gluon propagator 
as 



p6 + p4 (^^2 + 2M2) + p2 (2m2M2 + M4 + - ppt) + (M^ - ppt) + M^A^ - f (p + pt) 

(2) 

where the condensates are given by 

{A;A';)^-m' (^M)^P (3) 

and A'^ is related to the Gribov parameter 7 through A^ = 2g'^N'y^. We note that, when 
{(Pi^Pi) = {^i'<^i) then p = p"!", and the propagator reduces to 

yP ^ p4 + p2 (^2 + ^2 + + ^2 (J\^2 + p^) + A4 ^ ^ 

These propagators have the nice property to have -D(O) 7^ and so, they agree with lattice 
data. One of the most interesting properties of these propagators is that they can be 
rewritten as a sums of Yukawa propagators. So, for eq.([2]) one has 

Dip') = + + ■ (5) 

and for eq.dl]) is 

^ P = + -T—^ ■ 6 

This kind of propagators are in agreement with the scenario recently emerged from lattice 
computations. We will discuss this in the following sections showing how all this agrees with 
the current understanding of the Yang- Mills scenario on the lattice. 



B. Gluon propagator at the trivial infrared fixed point 

Computations on the lattice for a pure Yang-Mills theory showed that the running cou- 
pling manifests a peculiar behavior: The infrared limit is seen to reach zero lowering mo- 
menta 26 1 . 
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FIG. 1: Running coupling for 64^ and 80^ at /? = 5.7 taken from Ref. 



26| 



(http://www.sciencedirect.coin/science/jouriial/03702693 courtesy of Andre Sternbeck). 



From figiU it is easily realized that the coupling reaches trivial fixed points at higher and 
lower energies reaching a maximum in the intermediate regime. It must be emphasized that 
this is only true for a pure Yang-Mills theory. In QCD, the presence of quarks makes this 
fixed point non trivial. A trivial fixed point implies a Gaussian leading order functional for 
the partition function and we could write 



This kind of generating functional can be understood if we assume, as done in '80 [27[, that 
at low energies the functional of a Yang-Mills theory admits a series expansion in powers 
of the currents. Lattice data confirms a posteriori this hypothesis. Now we will show that 
exists a set of solutions for the inhomogeneous equations of motion of the Yang-Mills theory 
that admits a propagator of a free theory. Before to accomplish this, we note that, f or g auge 



theories, one can devise a reformulation of Wightman axioms as done by Strocchi 28[ that 



provides a Kallen-Lehman representation with a non-positive spectral function. If this is 
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the case, the most general form of propagator for a free theory takes the from 



n=0 



p2 — + ie 



(8) 



being m„ an internal spectrum of the free particle. Indeed, to put in Gaussian form the 
generating functional, we have to show that exists at least a gauge choice for classical Yang- 
Mills equations having as a solution eq.(IH]). So, following a standard notation (29|], we take 
for the generators 

= z/^'-^t- (9) 

being /"^'^ the structure constants of the gauge group. Then, for the covariant derivative we 
have 

D, = d,- igf^Al (10) 
with g the coupling constant and A° the potential. This gives 



(11) 



So, the equation of motion are 



(12) 



being arbitrary currents. To these equations we can add a gauge fixing term dy{d ■ A°')/S,- 
Now, we can prove that exist a solution of these equations with a propagator given by eq.(IHl). 
Taking into account eq.l ITT]) one has 

d^d^Al-(^ - d,id-A'')+gr''^d'^iAlA:)+gr''^A'^'id^Al-d,A';^)+g^^ = - 

13) 



Our aim is to look for a set of instanton solutions 



30| that can be cast into the form 31 1 



K = vim 



(14) 



being ?7jJJ some coefficients that can carry some dependence on momenta. Then, by a direct 
substitution into eq. (fT3!) we get 

1 
1 



1 - - 1 9.(r/'^-a0)+^7/'^SJr/ea'^(02)+^7/'^S>(r/e9^0-r/^9.0)+//'^^7'^''^r/'^S;^:^ 



(15) 
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By multiplying by rj""' both sides we arrive at 

V''''Vud'd,<P - (^1 - r/'^ . 9(r/'^ ■ d<j,) + g'v'"'r'''r''v'Xvt<P' = -J- (16) 

Now, we note that t]'^""!]'^ = N"^ — 1 (one can see this immediately by considering the SU(2) 
case) and the equation reduces to 

d''d,<P - [l - (r/-^ ■ dfct^ + = -j^. (17) 

where we have set = j /{N"^ — 1). Reducing to the case of the Lorenz gauge (^ = 1) being 
this equivalent to the Landau gauge at the classical level we are working on, this equation 
becomes classically exact (otherwise our solutions will be just asymptotic ones 2J]). The 



solution of this equation has been extensively discussed ( 32j and Refs. therein) and admits 
an iterative solution having the form at the leading order 



j^,m.-,)ui^) (18) 



with 



D{P') = E 2 2^- (19) 



n=0 

and 

= (2" +l)7p(i) M, ■ (20) 

being K{i) = j^^ ^ '^^^ ^ ^ 1.3111028777, and a formula for the spectrum of the theory, in 
the strong coupling limit, given by 

m„ = (2n + 1) ( ^1 ' A. (21) 

From this "mass spectrum" we can identify a string tension when we set 

^ = (^) ' A = {2rrNas)-^A. (22) 
Here A is just an arbitrary parameter arising from the integration of the equations of motion 



of the classical theory. It is not difficult to see, using Callan-Symanzik equation 33[, that 
this propagator implies P{Ng^) = ANg"^ and so the theory is infrared free, making this 
argument completely consistent with a Gaussian generating functional. 

So, there exists a class of classical solutions for the Yang-Mills equation in the Landau 
gauge that are consistent with the lattice results of a trivial infrared fixed point. We note the 



essential point that the sum of Yukawa propagators we obtain in eq. (fT9|) are exponentially 
damped and so, a few terms are already enough to fit a lattice propagator in the deep 
infrared limit. Our aim will be to compare it with a RGZ scenario and show that we are 
describing the same physics. This will provide fixed values to the free parameters, being 
these the condensates, of the RGZ scenario. 



III. COMPARISON WITH RGZ SCENARIO 

From eq. (l20|) we can immediately evaluate the weight of each propagator into the series 
([in]). It is easy to see that 

Zo^l.U Zi ^ -0.15 Zaf^O.Ol Z3 ^ -6.7 ■ 10'^ . . . (23) 

and so, taking the first few terms is enough to get a proper approximation in the deep 
infrared regime. Firstly, we fix the parameters for the case ([2]). We will have from eq. fll9p 

+ p^(mf Zo + mlZo + mlZi + mlZi + m^Z2 + rnfZ2) + mfmg^o + mlmlZi + ■m^m\Z2 



(24) 

where we have extracted an overall constant C = + + Z2 ~ 1. This gives, by a direct 
comparison with eq.(l2]), 

2M^ = m\ZQ + m\ZQ + m^Zi + m\Zi + m'^Z2 + m\Z2 
— pp^ = rn^m^ZQ + rn^m\Zi + rn^m\Z2 



+ 2M^ = ml + ml + ml 



2m^M^ + + - pp^ = mlml + mlmj + mjml 
m^ (M^ - ppt) + M^X^ - ^ (p + pt) = mlmfml (25) 



we 



In order to evaluate these parameters we have to fix the masses m„. As chosen in Ref. 
take ^/o^ = 0.44 GeV for the string tension. This will give, using eq.f l2ip . m^ = 0.527 GeV 
mi = 1.58 GeV and ma = 2.63 GeV. We get 



M ^ 2.2 GeV pp^ ^ 4 GeV^ = 0.0033 GeV^ (26) 



and, already at this stage, we note that, using the following formula for the condensate 



(g'A^) = -j^^^' (27) 
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one should have < 0. So, for consistency reasons we have to consider just the first couple 
of Yukawa propagators in eq. f|T9l) giving 



with 



2 

D{p') = {Zo + Z,) ^ f (28) 



S = M^+p,= ^0^1 + ^^^0 ^2^^^^2^^2^^2 ^2^^2^^;^4^^2^2 ^39) 

and in this case we get 

Zo + Zi^ 0.99 s ^ 2.85 GeV^ u ^ 1.67 Gel/ t ^ 0.83 Gel/^ (30) 

that gives 

{g'A') = -^^^rn' - 0.13 GeV'. (31) 
This gives a substantial validation to RGZ scenario for a gluon propagator in the form 



of eq.(!4]). Similarly, our scenario, as depicted in [32| is perfectly consistent with RGZ 
scenario and provides further validation with respect to lattice data. Similarly, we have 
A4 = 0.89 GeV^. Finally, D{0) = M^/\^ = 5.44 GeV~^ that is of the right magnitude as 
expected from lattice computations. 

We note that, while the value of the other parameters is perfectly consistent with other 

n □ 

determinations given in |l6l . llTj . the value of the dimension two condensate, {g'^A^), is 
significantly smaller. The reason for this can be traced back to the dependence on the 
coupling and on the lattice spacing. In fact, the determination in JigI . IitI is obtained from 
lattice data. Fitting at /3 = 6 and /3 = 2.2 provides quite different values for this condensate 
while we note that we are working at the trivial infrared fixed point however confirmed by 
RGZ scenario. This can be seen in a quite straightforward way. Using eq.( l2T]) for the mass 
spectrum we have 

2 2,2 2,2 Zprnj + Ziml \Zo + 9Zi\ TT^ r- 2 /ooN 

m = nin + m, — s = rrir, + m-, = 1= VNgA (32) 

° ' ° ' Zo + Zi Zo + Zi 4v^A'2(z) ^ ^ ^ 

and it is < being Zi < and 9\Zi\ > Zq. So, we note that, for lattice computations, 
are critical both coupling and lattice spacing and this explains the observed discrepancy in 
{igI . Iirj l . Similarly, such a conclusion can also be drawn for + pi and and we have 



Zq + Zi 4V2K^{i) 



Zo + Zi \ Z(j + Zi 



^ --,2 A 4 



(33) 



-Ng'A' 



This gives a proof of existence from first principles of the condensates in the RGZ scenario 
and, as an aside, we were able to fix them. 

IV. CONCLUSIONS 

In this paper we were able to show how RGZ scenario is perfectly consistent with a 
set of instanton solutions of Yang- Mills equations in the Landau gauge. The values of the 
condensates are finite and non-null making all the scenario well sound and explaining the 
excellent agreement with data on lattice. As an aside, we were able to fix in a closed form 
the condensates. 

It is interesting to note that one of the first principles we started from to derive the 
condensates was the emerging data from lattice proving the existence of a trivial infrared 
fixed point for a pure Yang-Mills theory. This gives a Gaussian generating functional with 
the propagator having all the parameters fixed by the theory. With a proper set of Wightman 
axioms, it is possible to get a proper extension of the Kallen-Lehman representation also for 
a gauge theory and this proves that the propagator one gets from the instanton solutions we 
obtain is indeed the one of a free theory. RGZ scenario also shows that the gluon propagator 
is the sum of Yukawa propagators and so both views validate each other. 

As a final consideration, it is interesting to note that, with all the parameter properly 
fixed, RGZ scenario can be used to work out a lot of phenomenology. A gluonic Yukawa 
potential is enough to produce a Nambu-Jona-Lasinio model from QCD lagrangian. 
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